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Abstract 

Sampling theory concerns the problem of reconstruction of functions 
from the knowledge of their values at some discrete set of points. In 
this paper we derive an orthogonal sampling theory and associated 
Lagrange interpolation formulae from a family of bounded rank-one 
perturbations of a self-adjoint operator that has only discrete spectrum 
of multiplicity one. 
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1 Introduction 



Sampling theory is concerned with the problem of reconstruction of functions, in a point- 
wise manner, from the knowledge of their values at a prescribed discrete set of points. 
Resolution of concrete situations in this theory generally involves the characterization of 
a class (usually a linear set) of functions to be interpolated, the specification of a set of 
sampling points to be used for all the functions in the given class, and the derivation of an 
interpolation formula. 

The cornerstone for many works on sampling theory is the Kramer sampling theorem 
[TO] and its analytic extension [6] . Orthogonal sampling formulae often arise as realizations 
of this theorem. The celebrated Whittaker-Shannon-Kotel'nikov sampling theorem [9J [171 
22J is also a particular case of the Kramer theorem, although historically the former came 
first and motivated the latter. 

Orthogonal sampling formulae have been obtained in connection with differential and 
difference self-adjoint boundary value problems (see for instance [21] and, of course, the 
paper due to Kramer himself [10]), and also by resorting to Green's functions methods 
[U [23], among other ODE's techniques. These results suggest that the spectral theory 
of operators should provide a unifying approach to sampling theory. Following this idea, 
a general method for obtaining analytic, orthogonal sampling formulae has been derived 
in [15] on the basis of the theory of representation of simple symmetric operators due 
to M. G. Krein [UJ [121 131 EE]- Roughly speaking, the technique given in [15] consists 
in the following: By [TT1 [T2l [T5] , every closed simple symmetric operator A in a Hilbert 
space 7i generates a bijective isomorphism between 7i and a space of functions H, with 
certain analytic properties. If the operator A satisfies some additional conditions, all of its 
self-adjoint extensions have discrete spectrum and every function / in Ti can be uniquely 
reconstructed, as long as one knows the value of / at the spectrum of any self-adjoint 
extension of A. 

Loosely speaking, the self-adjoint extensions of a symmetric operator with deficiency 
indices (1,1) constitute a family of singular rank-one perturbations of one of these self- 
adjoint extensions [31 Sec. 1.1-1.3]. Therefore, the methods developed in [15] also holds 
for singular rank-one perturbations, provided that these operators correspond to a family 
of self-adjoint extensions of some simple symmetric (hence densely defined) operator. We 
cannot use, however, a family of bounded rank-one perturbations in applications to sam- 
pling theory without making substantial changes to the results of [15]. Krein's approach 
to symmetric operators with equal deficiency indices does not work for bounded rank-one 
perturbations since operators of this kind may only be seen as self-adjoint extensions of a 
certain not densely defined, Hermitian operator [31 Sec. 1.1]. 

The main motivation of the present work is to develop a method in sampling theory 
analogous to [J5] for the case of bounded rank-one perturbations. With this purpose in 
mind, we begin by constructing a representation space for bounded rank-one perturba- 
tions in some sense similar to that of Krein for simple symmetric operators. Elements of 
this representation space are the functions to be interpolated. Then we obtain a general 
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Kramer-type analytic sampling formula which turns out to be a Lagrange interpolation 
formula. We also characterize the space of interpolated functions as a space of meromor- 
phic functions with some properties resembling those of a de Brange space. Examples are 
discussed in the last part of this work. 



2 Preliminaries 



The following review is based on the standard treatment concerning rank-one perturbations 
of self-adjoint operators, as discussed in detail by Gesztesy and Simon [HI [18], and Albeverio 
and Kurasov [3J. 

In a separable Hilbert space H, we consider a possibly unbounded, self-adjoint operator 
A with discrete spectrum of multiplicity one. Let \i be a cyclic vector for A, that is, 
{(A — ziy 1 ^ : z G C} is a total set in 7i. Throughout this work we assume that = 1. 

Given /i, let us define the family of bounded rank-one perturbations of A, 

A h -.= A + h(fi r )fji, heR, (2.1) 

where the inner product is taken, from now on, anti-linear in its first argument. Naturally, 
Dom(Ah) = Dom(A) for any h G R. Elementary perturbation theory implies that all the 
elements of (12 .11) have discrete spectrum. As pointed out in [3j Sec. 1.1], the operators 
Ah may be seen as self-adjoint extensions of some Hermitian operator, in the sense of |21j . 
with non-dense domain. 

Consider the family of functions 

F h {z) := (//, {A h - ziy 1 ^) , z <£ Sp(A) , heR. (2.2) 

In the sequel we shall denote F by F. By the spectral theorem, Fh(z) is the Borel transform 
of the spectral function mh(t) = (//, E h (t)fi), where E h (t) is the spectral resolution of the 
identity corresponding to Ah- Hence, Fh(z) is a Herglotz meromorphic function having 
simple poles at the eigenvalues of Ah. 

From the second resolvent identity [2TI Thm. 5.13] one obtains 

(Ah-ziy 1 = (a - ziy 1 - h((A h - ziy 1 ^, -}(a - ziy 1 , heR. 

This equation yields the well-known Aronzajn-Krein formula [181 Eq. 1.3] 

The spectral properties of the whole family (12. ip are contained in (12. 3p . Indeed, one can 
easily show that the function F h /F h i, h ^ h', is also Herglotz and its zeros and poles are 
given by the poles of Fy and Fh respectively [16]. Thus, the spectra of any two different 
elements of the family (12 .ip interlace, i.e., between two neighboring eigenvalues of one 
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operator there is one and only one eigenvalue of any other. Also, (12.31) implies that Xq G R 
is a pole of F h if and only if 

jk + h ^ (2 ' 4) 

Therefore for any i6l which is not a zero of F, there exists a unique h G R such that 
x is an eigenvalue of A^. One can extend this result to every x G R by considering an 
infinite coupling constant h = oo in (12.11) (see [HI Sec 1.5], [31 Sec. 1.1.2]). From the 
properties of F(z), it is shown that A ro also have simple discrete spectrum (see footnote in 
[U p. 55]) and Sp(Aoo) = {x G R : F(x) = 0}. Thus, for any x G R, there exists a unique 
/i£lU {oo} such that x is an eigenvalue of A h . 

The main peculiarity of A^, that separates it from the family A^ with finite h, is its 
domain. Indeed, for the kind of perturbations considered here, the domain of A^ is the 
set {ip G Dom(A) : (<p,fj) = 0} Sec. 1.1.1], [18, Thm. 1.15]. 

3 Sampling theory 

Based on the theory of rank-one perturbations, we construct in this section a linear space 
of meromorphic functions and derive an interpolation formula valid for all the elements 
in H^. 

In our considerations below, the following vector-valued function of complex argument 
will play an important role. 

Notice that £(z) is well defined for z G Sp(A). 

Lemma 1. For any x G R \ Sp(Aoo) ; there exists (a unique) h G R suc/i i/wi 

C(i) e Ker(A ft - a;/) . (3.1) 

Similarly, 

(A - xrpV e Ker(A DO - x/) (3.2) 

/or every x G Sp(Aoo). 

Proof. We first consider x G" Sp(A) U Sp(A 00 ). Let h ^ be such that x G Sp(/U) (we 
already know that there is always such h). We have 

*««) = - «o- v = + ft) /. + ^ - «o-v . 

The first assertion of the lemma follows from the last expression and ( 12.41) . When x G 
Sp(A), the statement follows by a limiting argument based on the fact that A is closed. 
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We now prove the last assertion of the lemma. Define P := (//, •)//. By virtue of 
Thm. 1.18, Rem. 2], there is a cyclic vector 77 that obeys 

- ziyS = -^-(1 - P)(A - Z )-V , 2 £ s P (A) u Sp(A») • (3.3) 

Clearly, 77 G Dom(A 00 ). We compute the projection of 77 along the eigenspace associated 
to x G Sp(A 00 ). Using (13.31) and some e > sufficiently small, we obtain 

[E oo (x + 0)-E oo (x-0)]r l = ^- [ {A^-zY^dz 

2m J\x-z\=e 



\x — z\=e 

p, dz 



\x—z\ = € 



F{z 



1 (A-zI)- 1 -! 



Res — ]— ) (A - xl) V . 



z=x 



□ 



F(z) 

Since the last expression is different from zero, (13.21) is proven. 

Definition 1. For any ip G 7i, let $ M be the mapping given by 

:=<£(*),¥>>, «eC\Sp(A»)- 

We sometimes shall denote by 

The mapping $ M is a linear injective operator from H. onto a certain space of meromor- 
phic functions TC^ := $ M 7i. The injectivity may be verified with the aid of (13.11) . Some 
properties that characterize the set Ti^ will be accounted for in the next section. 

Proposition 1. Given some fixed fee R, let {xj}j = Sp(Af l ). Define Gh(z) := l/F^z) = 
h + 1/F(z). Then, for every f(z) G Ti^, we have 

/(*)= E (z JtfG'(x) f{Xjh ^^\MAoo). (3.4) 

The series is uniformly convergent on every compact subset of the domain. 

Proof. Because of the assertion (13. ip of Lemma [fl is a complete orthogonal set in 

TC. Hence 

<p(z) = m,<p)= E ( 3 - 5 ) 

where the series converges uniformly on compacts of C \ Sp(Aoo) by virtue of the Cauchy- 
Schwarz inequality. 
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Now, the first resolvent identity implies 



(£(zUH) = (s-BJ) 



_F(w) F(z) 



In conjunction with (12.41) and the convention l/oo = when h = 0, the last equation gives 
rise to the identity 



F(*) 



Finally, notice that 
G' h (w) = 



1 



-F'(w) 



(lM, {A-V,)*VL) = -(t(lD),t@)). 



F 2 (w) v ' F 2 (w) 
Evaluation of the last expression at w — Xj yields the desired result. 



□ 



Remark 1. Equation ( 13. 5ft is an orthogonal sampling formula of Kramer-type (TO]- Since 
the function Gu{z) has simple zeroes at the points of Sp(v4/ l ), expression (13 .4p is indeed a 
Lagrange interpolation formula. 



4 Spaces of interpolated functions 

The present section is devoted to the characterization of the space of functions intro- 
duced by means of the mapping $ M of Definition [TJ Notice that TC^ depends on both the 
operator A and the cyclic vector /i. 

The following statement gives a quite explicit description of Ti^. 

Proposition 2. 



^ M = < f(z) = c+ Y : c, c n <£ C, Y \c n \ 2 F'(x n ) < oo 

where the series above converge uniformly on compact subsets ofC\ Sp(A 00 ). 

Proof. Let G denotes the set defined by the right-hand side of the statement. 

Given x n G Sp(A OQ ), it follows from (13. 2p of Lemma [T] that u(x n ) := (A — x n I)~ l n is 
the associated eigenvector (up to normalization). Taking into account the first resolvent 
identity, we get 



u(x n )(z) = -J--({A - z) 1 n,(A-x n ) V) 



F(z)" ' ^ " F(z) 



^ + F(x, 



Z Xr 



Z Xr 
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The expression after the first equality above also implies that ||a;(a; n )|| 2 = F'(x n ). Recalling 
the definition of Dom(A 00 ), it follows that the set 

B := {//}U {IK^irM^L^sp^) 

is an orthonormal basis in 7i. Now, take an arbitrary element ip of TC and expand it on 
the basis B. By applying $ M to ip and using the Cauchy-Schwarz inequality, we conclude 
that C G. 

The inclusion G C follows from noticing that any f(z) £ G is the image under $ M 
of an element in TC of the form cfi + SznsSpfAjo) c nu(x n ). 

□ 

Notice that, by virtue of Proposition [2} the only entire functions in TC^ are the constant 
functions. Also, one easily verifies that any constant function in Tt^ is the image under $ M 
of a vector in Span{/i}. 

The following straightforward result shows that the functions in TC^ share some prop- 
erties with those in a de Branges space. 

Lemma 2. The space TC^ has the following properties: 

(i) . Assume that f(z) £ TC^ has a non-real zero w. Then g(z) := jz^f( z ) a ^ so belongs to 

(ii) . The evaluation functional /(•) i— > f(z) is continuous for every z £ C \ Sp(A 00 ). 
(Hi). For every f(z) £ TC^, g(z) := f(~z) belongs to TC^. 

Proof. We have f(z) = (£(z),</?) for some ip £ TC. Given w such that f(w) = 0, consider 
i] = (A — wI)(A — wl)~ l ip. A short computation yields g(z) = (£,(z), rj), thus showing (i). 
Assertion (ii) is rather obvious so the proof is omitted. On the basis of Proposition [2] one 
verifies (iii). 

□ 

In what follows we show that TC^ can be endowed with several Hilbert space structures, 
each one determined by the spectral functions rrih(x), h £ R. 

Lemma 3. Let h £ R and {xj}j = Sp(Ah), arranged in non-decreasing order. Then the 
spectral function rrih(x) is given by 

Proof. Let us recall first the following well-known results [T8l Thm. 1.6] 

limeReF/ l (x + ie) = 0, 
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hmelmi^z + ie) = gm^ ^-^L = m/i({a;} ) , 
and also the identity 

(A - z/) ~ V = T77b^- z/) ~ 1 ' 1 ' 

Consider Xj G Sp(A^). It suffices to verify that rrih({xj}) = \\(,(xj)\\~ 2 ■ By resorting to 
the equalities mentioned above, a straightforward computation shows that 

(£{xj - ie),£(xj - ie)) = 1 2 ((A - (xj + ie)/)"V, - (a;, + ie)I)~V) 



\F h ( Xj + ie)\ 2 



e 2 dm h (y) 

[eReFhixj + ie)f + [elmF^x, + ie)} 2 J R (y - Xj ) 2 + e 2 
1 



e -> 0. 



The proof is now complete. 

□ 

Remark 2. For h ^ 0, this result is in fact statement (ii) of jTHJ, Thm. 2.2] in disguise. 

Proposition 3. For arbitrary feel, the map 3^ zs a unitary transformation from 7i onto 
L 2 (R,dm h ). 

Proof. $ M is a linear isometry from 7i into L 2 (IR, drrih)- Indeed, 



life) II 



xje8p(A h ) 

Now consider f(x) G L 2 (IR, dm^). This means that 



1= E M§<- 



Sj6Sp(yl h ) 

Define 



iifeoir 



xjeSp(A) 



llfe)ll' 
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which is clearly an element in Ti. It is not difficult to verify that ||/(-) — ^( - )IL = 0- 

□ 

Define C := fc" 1 ^, where (C 'f)(z) = 7M for / G H„. By (iii) of Lemma [2] and 
Proposition [3l it follows that C is a complex conjugation with respect to which both A 
and \i are real. 

We conclude this section with a comment about the representation of the operators 
as operators on Ti^. A simple computation shows that, for every h eR, A h is transformed 
by $^ into a quasi-multiplication operator, in the sense that 

Aw{z) = -=- r r(fJi,<p) + z<p(z) (4.1) 

for every if G Dom(A). This is obviously the multiplication operator in <3> M Dom(yl 00 ). 
Moreover, (14.11) reduces to the multiplication operator in a weak sense; indeed, 

&(-)) h = (?(.),(.)$(■)) h , 

for every if, ip G Dom(y4). 



5 Examples 



Rank-one perturbations of a Jacobi matrix. 

Jacobi matrix 





bi 








h 


1i 


b 2 








b 2 


?3 


63 












V- 









Consider the following semi-infinite 



(5.1) 
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with q n G R and b n > for n G N, and define in the Hilbert space / 2 (N) the operator J 
in such a way that its matrix representation with respect to the canonical basis {5 n }^ =1 
in l 2 (N) is (15.11) . By this definition, J (cf. [2j Sec. 47]) is the minimal closed symmetric 
operator satisfying 



(5 n , J5 n ) = q n , (5 n+ i, J5 n ) = (S n , J5 n+ i) = b n 



Vn G N. 



The Jacobi operator J may have deficiency indices (1,1) or (0,0) [1, Chap. 4 Sec. 1.2], 
[T91 Cor. 2.9]. For this example we consider J to be self-adjoint, i. e., the case of deficiency 
indices (0, 0). We also assume that J has only discrete spectrum. Our family of self-adjoint 
operators is given by 

J h := J + h(5 1 ,-)5 1 , heR. (5.2) 
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It is relevant to note that 5i is a cyclic vector for J since the matrix elements b n are always 
assumed to be different from zero. 

One can study J through the following second order difference system 

bn-lfn-l + Qnfn + Kfn+1 = zf n 71 > 1 , Z E C . (5.3) 

with boundary condition 

+ bih = zfi ■ (5.4) 

If one sets fi = 1, then f 2 is completely determined by (j5.4p . Having fi and f 2 , the 
equation (15. 3p gives all the other elements of a sequence {f n }^=i that formally satisfies 
(15.31) and (15. 4p . f n is a polynomial of z of degree n — 1, so we denote /„ =: P n _i(z). 
The polynomials P n (z), n = 0, 1, 2, . . . are referred to as the polynomials of the first kind 
associated with the matrix (15. ip . The polynomials of the second kind Q n {z), n = 0, 1, 2, . . . 
associated with (15. ip are defined as the solutions of 

K-lfn-l + Qnfn + Kfn+1 = Zfn TL E N \ {1} 

under the assumption that fi = and f 2 = b^ 1 . Then 

Qn-l(z) :=/n, VnGN. 

Qn( 2 ) is a polynomial of degree n — 1. 

Let -P(z) = {P n (^)}^L and Q(z) = {Q„(z)}^L . Then, classical results in the theory 
of Jacobi matrices [Tj give us the following expression for defined in Section [3J 

= ^(*) + ^G(z) > 

where F(^) is the function given by (12.21) with h = 0. In this context P(z) is referred to 
as the Weyl function of J and may be determined by 

F(z) = - lim — ^ , w n (z) := £f\ (5.5) 

where the convergence is uniform on any compact subset of C \ Sp(J) pfl Sees. 2.4,4.2]. 

The operator Joo corresponds in this case to the operator in Z 2 (2, 00) whose matrix 
representation is (15.11) with the first column and row removed. 

For any / E TCs 1 there is a sequence {(Pk}kLi e ^ 2 (N) such that 

f(z) = [Pk-i(z)<p k + -^rQk-i(z) 
fc=i v r w 

Notice that the poles of / are the eigenvalues of J^. 
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By Proposition [T] we have the following interpolation formula 

/(*)=lim £ h ~ w n( z ) /( ) heR . (5.6) 

Indeed, one can write (15. 6p on the basis of (13. 4p using the uniform convergence of the 
limit and the series in (15.51) and (13.41) . respectively, and the fact that w' n {z) is also uniform 
convergent. 

One-dimensional harmonic oscillator. In this example we look at the sampling formula 
provided by only the unperturbed operator. 

On L 2 (R, dx), consider the differential operator 

A ^2 

A := —3—^ + x , 

dx z 

which is essentially self-adjoint on C^°(R). The eigenvalues are 2n + 1 for neNU {0}; the 
corresponding eigenf unctions are 

<f> n (x) = n-^(2 n n\)-^ 2 e~ x2 / 2 H n (x), 

where H n (x) are the Hermite polynomials. 
A cyclic vector for the operator A is 

00 

"M = E wW«*> = ^e-i<^^) , (5.7) 

n=0 V 

where the last equality follows from a quick look to the generating function of the Hermite 
polynomials. 

Now, for every z Sp(-A), 

00 1 

m = £ 



n\(2n + 1 - z) ' 

n=0 V ; 

An elementary argument involving series of partial fractions then shows that 



F(z) = -/ e- cos0+lsm0 e l — e d9. (5.8) 

2 cos — J_ n 



(See, for instance, [5].) 

The mapping $ M is defined by the vector-valued function 



[(A-^)-y](x) 1 r° 

F(z) F(z) 
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where the Green's function K(z;x,y) is given by (see [20] ) 



K(z;x,y) 




TV 



1/2 



cosf 



X 



%(2 



D,-i(- 



1/2 x)D,_ 1 (-2 1 / 2 y), y<x, 
-2 1 / 2 x)D^{2 1 / 2 y) , y>x. 



(5.9) 



In the last expression, D p (x) denotes the parabolic cylinder function of order p. 
In a fashion more customary for sampling theory, we state the following result: 

Corollary 1. Let /i(x), F(z) and K(z; x,y) be given by ( |5. 7| ], 115. 8\) and Ii5.9\) . respectively. 
Let G(z) := 1/F(z). Then, every function f(z) of the form 
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